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Light propagation in coupled optical waveguides provides an excellent laboratory 
tool to visualize with classical waves the dynamical aspects embodied in a wide va- 
riety of coherent phenomena encountered in atomic, molecular or condensed-matter 
physics pQ. Among others, we mention adiabatic stabilization of atoms in ultrastrong 
laser fields [2], Bloch oscillations, Zener tunneling and dynamic localization in crys- 
talline potentials [2H5], coherent population transfer in laser-driven atoms [6j[7], and 
quantum Zeno dynamics [8]. As compared to quantum systems, optics offers unique 
possibilities such as a direct visualization in space of typical ultrafast phenomena in 



time and the ability to explore coherent dynamical regimes of difficult access in quan- 
tum systems or obscured by detrimental effects such as dephasing and particle-particle 
interaction. Previous optical analogies [1] deal with coupled channel waveguides in an 
assigned geometry, in which spatial light transport turns out to be analogous to the 
coherent evolution of a single-particle electronic wave function. Unfortunately, these 
photonic structures could not account for vibrational effects typical of e.g. molecules or 
crystals. Nuclear degrees of freedom in molecules are known to yield a wealth of inter- 
esting phenomena, especially when driven by ultrastrong laser fields. Coherent effects 
in molecules, like laser-assisted preparation and control of vibrational wavepackets 
and molecular bond engineering, provide theoreticians and experimentalists with ex- 
amples of time-dependent quantum mechanics at work, with important applications 
to femtochemistry [9HTT]. Though current pump-probe techniques enable the obser- 
vation of time-resolved wave packet dynamics of nuclear motion, optical analogues 
of coherent vibrational phenomena may offer the possibility to assess wave packet 
dynamical phenomena of difficult access in the quantum realm, such as wave packet 
dynamics and momentum filtering at surface potential crossings [T2] . 
In this Letter optical analogues of laser-assisted coherent vibrational phenomena in 
molecules, such as light-induced molecular stabilization and wave packet dynamics 
at surface potential crossings, are theoretically proposed for light beams propagating 
in coupled slab waveguides. The photonic structure, schematically shown in Fig. 1(a), 
consists of two slab waveguides Si and S2, placed at a distance d, which confine the 
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field in the transverse y direction. In the weak guidance approximation, light prop- 
agation at wavelength A is described by the Schrodinger-like wave equation for the 
electric field envelope ip(x,y,z) [TIIT3"] 

%W z i) = -[h 2 /(2n s )}V 2 t ij + [n s -n(x,y)}ij (1) 

where V 2 = d 2 + d 2 is the transverse Laplacian, Ti = \/{2n) = 1/k is the reduced 
wavelength, n s is the bulk refractive index, n = n s + Arii(x,y — d/2) + An 2 (x,y + d/2) 
is the refractive index profile of the structure, and Ani^^, y) are the refractive index 
increase in the two slab guiding regions. The index profiles An\ y2 are generally allowed 
to slowly vary along the horizontal x direction, such as in case of graded-index slabs or 
slabs with slowly-varying thickness [13]. We indicate by qi(y; x) and g 2 (z/; %) the funda- 
mental mode profiles of the two slabs at a given transverse coordinate x and by Wi >2 (x) 
the respective effective indices, i.e. — (h 2 /2n s )d 2 qi i2 + Arii 2 (x, y)q\ y2 = Wi,29i, 2- As- 
suming that Wi t 2 {%) differ from a reference value Wq by a small quantity (of order 
~ e 2 ) and that the field ip and the effective indices W\ j2 var y slowly with respect to x, 
a multiple-scale asymptotic analysis shows that the solution to Eq.(l) can be written 
as ip(x, y, z) = [ipi(x, z)q 1 (y - d/2; x) + ifa(x, z)q 2 (y + d/2; x)} exp(-iW z/h) + 0(e 2 ), 
where the slowly-varying amplitudes ^1,2 satisfy the coupled Schrodinger-like equa- 
tions 

ihd e i/>x,2 = -[h 2 /(2n s )}dl^ 2 + Wi, 2 (x)^i,2 + «(x)^ 2 ,i. (2) 

In Eqs.(2), we have set U 1>2 = W lj2 {x) - W + f dyAn 2)l q\ 2 + (h 2 /2n s ) J dy(d x q h2 ) 2 ~ 
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Wip (x) — Wo, k = f dyAniqiq 2 = J dyAn 2 qiq2 and the normalization conditions 
J dyq\ 2 = 1 have been assumed. In their present form, Eqs. (2) are analogous to the 
quantum mechanical equations describing the vibrational dynamics of a molecule onto 
two energy surfaces coupled by an external laser field (interaction strength k) in 
the Born-Oppenheimer and rotating-wave approximations (see, for instance, [9|[T2]). 
Note that the temporal variable in the quantum problem is replaced here by the 
spatial propagation distance z, the reduced nuclear mass by n s , the Planck constant 
by the wavelength A, and the inter- nuclear distance by the coordinate x. This anal- 
ogy enables to investigate in optics the classical analogues of a wide variety of wave 
packet phenomena encountered in molecular systems. We will limit here to discuss 
two of such analogies. The first one, shown in Fig. 1(b), is the phenomenon of laser- 
induced stabilization of a dissociating molecule, which is usually explained by the 
formation of light-induced adiabatic surface potentials U±(x), defined by the equa- 
tion U± - (Ui +U 2 )U± - k 2 +U1U2 = (see, for instance, P H4"IIT5] ) . Let us consider 
a diatomic molecule (like Hj) characterized by a bonding {Ui) and an antibonding 
(U2) potential curve. A dissociating vibrational wave packet is initially created from 
the ground vibrational state at the internuclear distance X\ by a strong and ultra- 
short laser pulse at frequency U\ [Fig. 1(b)]. If a second continuous- wave laser field 
at frequency U2 is applied to set in resonance the two potential curves at a distance 
x 2 > Xi, a potential well is formed on the adiabatic surface U+(x) at the avoided 
crossing [see Fig. 1(b)], and molecular dissociation in thus inhibited in the adiabatic 
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limit. The nuclear distance x 2 at which the new bond is formed can be tuned by vary- 
ing the frequency u 2 of second laser. The optical realization of such a process is shown 
in Fig. 2 (a). Here the two slab waveguides have a slowly- varying thickness, linearly de- 
creasing with x for slab Si and linearly increasing for slab S2. The slabs are truncated 
at x = xq, providing a steep potential barrier for light. An equal slab thickness is 
reached at x = x 2 > xq. Correspondingly, the potential curves IA\ and U 2 entering in 
Eqs.(2) mimic the bonding and anti-bonding potential curves of Fig. 1(b), with a level 
crossing at x = x 2 . Moreover, the coupling rate k turns out to be nearly independent 
of x. In our waveguide system, excitation of a dissociating wave packet at x — x\ on 
the antibonding surface is simply achieved by launching an elliptic Gaussian-shaped 
beam into the slab S2 at normal incidence [inset in Fig. 2(a)]. Inhibition of molecular 
dissociation, arising from the formation of an adiabatic trapping well as the coupling 
strength k increases, is shown in Figs.2(b)-(d). The integrated intensity light dis- 
tributions J g \i/)(x, y, z)\ 2 dy oc \ipi(x,z)\ 2 and J" g \i/)(x, y, z)\ 2 dy oc \ip 2 {x 1 z)\ 2 trapped 
in the two slabs, shown in Fig.2(b)-(d), are numerically computed by a standard beam 
propagation analysis of Eq.(l). Note that, for a large waveguide separation [Fig.2(b)], 
the wave packet on the antibonding curve is accelerated toward increasing values of x 
and the optical analogue of molecular dissociation is attained. As the waveguide sep- 
aration decreases (and thus the coupling k increases), wave packet trapping is clearly 
observed [Fig.2(d)]. In this case, at the crossing point x = x 2 the wave packet period- 
ically tunnels between the two waveguides, and is accelerated in opposite directions 
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owing to the opposite slopes of the curves IA\ and U 2 - This trapping mechanism is 
imperfect due to Landau-Zener tunneling at each level crossing, and the lifetime of 
the trapped oscillatory motion increases as k increases [compare Figs. 2(c) and (d)]. 
The second analogy is the acceleration of a molecular wave packet after a potential 
crossing due to the dependence of the Landau-Zener transition probability on the 
wave packet velocity [12]. Let us consider a wave packet that propagates, at a con- 
stant speed t>o, along a flat potential curve U\{x) = and undergoing a Landau-Zener 
crossing with a second potential curve Uiix) = F(x — £2) at x = £2 [see Fig.l(c)]. 
In the semiclassical approximation, the wave packet behaves like a particle with a 
definite momentum p = n s Vo, which is assumed to weakly spread during the cross- 
ing and to keep its initial velocity vq. After the crossing, the probability Pi of the 
wave packet to remain on the original potential surface U\ is given by the usual 
Landau-Zener formula P[12] P\ = exp(— 7rA), where A = 2k 2 / {\F\Tivq) . Note that 
Pi is larger for a faster wave packet. If the wave packet momentum distribution is 
broad around its mean po, a more complex scenario is found. In particular, owing to 
the dependence of Pi on v , a filtering effect in the momentum space occurs, with 
the faster components of the wave packet preferentially remaining on the original 
potential curve. The final result is a slight increase of the mean wave packet veloc- 
ity [12] . The analogue of this phenomenon for light beams is shown in Fig. 3. Here 
two slab waveguides, the first one with constant and the second one with a linearly- 
increasing thickness [Fig.3(a)], mimic the energy crossing scenario of Fig. 1(c) with a 
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force F ~ 34 /im/cm 2 and a coupling k/Ti ~ 1.04 cm -1 . A tilted Gaussian-shaped 
wave packet i[)(x,y,0) = exp{ — [(y — d/2)/w y ] 2 — [(x — xi)/w x ] 2 — iv n s kx}, with 
w x = 8 /zm, w;^ = 2.4 /mi, <i = 10.7 /mi, and Xi = —800 /xm, excites waveguide Si 
with an initial transverse velocity (refraction angle) vq = 278 /im/cm, corresponding 
to a Landau-Zener probability Pi ~ 0.48. Figure 3(b) depicts the path followed by 
the center of mass (x) of ip\ (solid line); for comparison, the straight path followed 
by ipi in absence of slab S2 is also depicted (dotted curve). The two insets in Fig.3(b) 
show the the evolution of the fractional beam power P\{z) of wave packet ipi and 
its instantaneous velocity v(z) = d(x)/dz. Note that, as the wave packet reaches the 
crossing region and light transfer to S2 starts, the velocity v(z) first decreases from 
the initial value vq, then increases and reaches an asymptotic value slightly larger 
than v , indicating an acceleration of the wave packet. Such a behavior can be ex- 
plained [12] after observing that the faster wave packet components reach first the 
crossing region and partially tunnel into slab S2, so that a deceleration of ip\ is ini- 
tially observed. When the slower wave packet components reach the crossing regions, 
they are partially transferred into slab S2 and an acceleration of ipi is thus observed. 
Since the tunneling probability is smaller for the faster wave packet components, the 
final mean velocity of ip\ is larger than the launching one v . This effect, however, is 
small because the change of beam path remains internal to the diffraction cone of the 
beam [see the shaded area in Fig.3(b)]. Similar behavior is obtained by reversing the 
slope of the potential U2. 
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Optical analogs of coherent phenomena in matters have so far visualized the behavior 
of the electronic wave function solely pp. Here a photonic structure capable of mimick- 
ing coherent vibrational effects in molecules has been proposed. The present results 
may stimulate theoretical and experimental studies on optical analogs of coherent 
vibrational phenomena in molecules or crystals. 
Author E-mail address: longhi@fisi.polimi.it 
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List of Figure captions 

Fig.l (Color online) (a) Schematic of a two-slab optical waveguide structure, 
(b) Nuclear potential energy diagram of a diatomic molecule interacting with two 
laser fields, (c) Schematic of wave packet splitting at a linear potential energy crossing. 

Fig. 2 (Color online) (a) Refractive index of the slab waveguide structure that 
mimics the molecular energy diagram of Fig. 1(b). The inset shows the elliptical 
Gaussian beam that excites guide S2 at x = x±. (b)-(d): Evolution of light dis- 
tributions \ipi(x,z)\ 2 (right panels) and |^ 2 (x, ^)| 2 (left panels) for decreasing values 
of waveguide separation d: (b) d — 12 /mi, (c) d — 9 /mi, (d) d — 8 /mi. Other 
parameters are: A = 633 nm, n s = 1.45. 

Fig. 3 (Color online) Wave packet dynamics at a potential crossing, (a) Re- 
fractive index of the slab waveguides that mimics the potential crossing of Fig. 1(c). 
(b) Path followed by the center of mass of wave packet ip\ (solid curve). The shaded 
area shows beam spreading. The insets show the evolution of fractional beam power 
Pi (bottom) and speed v (top) of wave packet ipx. 
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Fig. 2. 



13 




14 



